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We show results from simulations of the Langevin dynamics of a two-dimensional scalar model with com- 
peting interactions for ultrathin magnetic films. We find a phase transition from a high temperature disordered 
phase to a low temperature phase with both translational and orientational orders. Both kinds of order emerge 
at the same temperature, probably due to the isotropy of the model Hamiltonian. In the low temperature phase 
orientational correlations show long range order while translational ones show only quasi-long-range order in a 
wide temperature range. The orientational correlation length and the associated susceptibility seem to diverge 
with power laws at the transition. While at zero temperature the system exhibits stripe long range order, as 
temperature grows we observe the proliferation of different kinds of topological defects that ultimately drive 
the system to the disordered phase. The magnetic structures observed are similar to experimental results on 
ultrathin ferromagnetic films. 

PACS numbers: 75.40.Mg, 75.40.Cx, 75.70.Kw 



I. INTRODUCTION 

In the last years the interest in understanding the thermody- 
namic and mechanical properties of magnetic ultrathin films 
has grown considerably. Part of this interest is obviously mo- 
tivated by the great amount of technological applications re- 
lated to their magnetic behavior, such as data storage and elec- 
tronics 1 . In recent years the advance in experimental tech- 
niques has made possible to study in great detail the complex- 
ity of magnetic order in thin films, where an extremely rich 
phenomenology is present 2 ^. Part of this phenomenology 
- as in many other physical and chemical systems 5 - is due 
to the presence of competing interactions acting on different 
length scales that frustrate the system and leads to mesoscopic 
pattern formation. In the case of an ultrathin magnetic film 
in the absence of an external magnetic field, stripe patterns 
of opposed magnetization are formed due to the competition 
between the ferromagnetic exchange interaction and the long- 
ranged dipolar interaction 6 -. 

Of particular interest for applications in data storage are 
films with strong perpendicular anisotropy, where spins point 
preferentially out of the plane. This happens when the (per- 
pendicular) anisotropy energy overcomes the effect of dipolar 
interactions which induce an in-plane anisotropy. When this 
happens the system goes through a "spin reorientation tran- 
sition" (SRT) 7 -^. In this phase, when magnetization points 
preferentially out of plane, complex magnetic structures arise, 
showing the formation of patterns with stripe orde r 2 ! 3 ' ° , Anti- 
ferromagnetic stripe order dominates the low temperature/low 
energy behavior. As temperature increases towards the SRT, 
perfect stripe order is disrupted by proliferation of topologi- 
cal defects, like dislocations and disclinations in the magnetic 
patterns, which eventually drive the system to the paramag- 
netic high temperature phase. Understanding how magnetic 
stripe order emerges from microscopic interactions and the 
characterization of the relevant thermodynamic behavior in 
the perpendicular phase is of relevance both from a funda- 
mental point of view and for future applications in magnetic 



memory devices at the nanoscopic level. 

Theoretical analysis of thermodynamic behavior rely on 
elastic Hamiltonian approximations and show a variety of pos- 
sible scenarios, strongly dependent on the behavior of dif- 
ferent kind of anisotropies present in the elastic energ y 6 ' 8,1 1 . 
Inli two possible scenarios were anticipated, in one there is a 
single phase transition from a high temperature paramagnetic 
phase to a smectic (stripe) phase at low temperatures. Other 
possible scenario shows the presence of an intermediate ne- 
matic phase, where translational order is lost but orientational 
one persists. Searching for evidence of these scenarios is one 
of our motivations for the present work. There is also a rather 
large number of numerical simulations and analysis of the dif- 
ferent patterns observed in ultrathin film models. Relevant to 
the present work are a series of simulations by Maclsaac and 
coworkers 12 . They made detailed Monte Carlo simulations 
of Heisenberg and Ising systems with competing exchange 
and dipolar interactions. The stripe nature of the ground 
state in a two dimensional system of Ising spins was deter- 
mined analytically in 13 . Phase diagrams showing the pres- 
ence of the spin reorientation transition in Heisenberg models 
with perpendicular anisotropy, dipolar interactions and with 
or without exchange interactions were also studied 9,14 . Can- 
nas and coworkers have made a series of detailed predictions 
on phase diagrams and dynamic properties of a two dimen- 
sional Ising model with competing exchange and dipolar in- 
teraction a 15 ' 16 ' 17 by Monte Carlo simulations. In particular, 
inii evidence was presented of a new, intermediate, nematic 
phase between the stripe and disordered phases. There is also 
a rather large amount of numerical work focused on qualita- 
tive descriptions of patterns, but few quantitative approaches. 

Simulations of Heisenberg models have been largely re- 
stricted to determinations of phase diagrams mainly because 
analyzing quantitatively the structure of phases and patterns 
that emerge is a computationally very hard task in these mod- 
els. On the other hand, studies of Ising models have also con- 
centrated on quantitatively determining phase diagrams and 
also the presence and evolution of magnetic patterns. Never- 
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theless in the Ising case the sharp nature of domain walls in- 
duces a rather artificial structure of patterns, where e.g. tetrag- 
onal symmetry dominates completely the scene in square lat- 
tices. In this work we introduce a model intended to make 
a bridge between the more realistic but cumbersome Heisen- 
berg models and the much simpler Ising models which nev- 
ertheless are not suitable for understanding the diversity of 
magnetic structures present in a real ultrathin film. Our re- 
sults are qualitative in the sense that we do not intend to repro- 
duce experimental parameter values of a particular system, but 
nevertheless we present a quantitative picture of phase transi- 
tions and the behavior of relevant thermodynamic variables in 
the regime of perpendicular magnetization, which can be eas- 
ily mapped to real situations. Furthermore, the model intro- 
duced below admits numerical as well as analytical treatment 
of both equilibrium and dynamical behaviors and a compre- 
hensive picture is beginning to emerge 1 ^. 

With this aim, in the present work we introduce and an- 
alyze, by means of Langevin simulations, a coarse grained 
model of an ultrathin film in two dimensions. By making rea- 
sonable approximations to the full micromagnetic dynamics, 
which should be satisfactory in the perpendicular region, we 
address the existence of antiferromagnetic stripe order, how 
thermal fluctuations affect this order, the possible appearance 
of a phase transition at finite temperature, and the characteri- 
zation of the magnetic structures relevant in each temperature 
region. Coarse-grained models similar to the one introduced 
in the next section have been studied in the past. In an early 
work, Roland and Desai 20 studied the dynamical process of 
phase separation and emergence of stripe order in Langevin 
simulations of a model for ultrathin films. More recently, 
Jagla 4 explored the different morphologies and patterns that 
can appear in such systems, showing a variety of very interest- 
ing phenomena. Furthermore, he showed that the same model 
is capable of reproducing the detailed phenomenology of his- 
teretic behavior known in ultrathin films 21 . Nevertheless, up 
to our knowledge, more quantitative studies of the phase tran- 
sitions and the different kinds of magnetic order present in 
these type of models have not been addressed up to now. 



n. MODEL AND NUMERICAL IMPLEMENTATION 

A widely acceptable microscopic description of micro- 
magnetic dynamics is given by the Landau-Lifshitz-Gilbert 



equation. 



— = —a 4> x B — 7 x {(f) x B) 



(1) 



where <ft is the three dimensional magnetic moment vector, 
B is the effective field acting on it and a and 7 are micro- 
scopic phenomenological constants. The first term induces a 
precessional movement of the magnetic moment around the 
field, while the second term is a phenomenological one repre- 
senting a damping effect which induces the magnetization to 
align with the field. This representation of magnetic dynam- 
ics is satisfactory in a wide variety of situations but it is very 



difficult to analyze analytically and is also very demanding 
computationally. Nevertheless, below the SRT where perpen- 
dicular anisotropy dominates, one can obtain a much simpler 
description of magnetization dynamics suitable to our pur- 
poses. Considering a single moment <\> with the effective field 
B pointing in the z direction and expanding the three compo- 
nents of (f> in eq.(fl]i, one easily finds that the evolution of the 
perpendicular component is given by 21 : 



dt 



1 B{l-<f>l) 



(2) 



where the constraint (\>\ + <fiy + <j) 2 z = 1 was used (this is auto- 
matically satisfied by the LLG dynamics). Then, in the limit 
of strong perpendicular anisotropy, when the effective field 
can be approximated to point preferentially along the z direc- 
tion, the evolution of the perpendicular component is approx- 
imately autonomous, although the other two components are 
not. Consequently from now on we will restrict the analysis 
to the z component and drop out the corresponding subindex. 

In the system of interest, the effective field will be com- 
posed typically of three contributions of the form: 



B = h + acj)-b 



J dV J(|x 



(3) 



where h is an external field, the second term is a perpendicu- 
lar anisotropy field and the last term will be the dipolar field in 
our case. The constants can be easily related with experimen- 
tal ones, but this is not necessary for the purposes of this work, 
which is to study some universal properties of the phases and 
patterns emerging in such a system. 

The complete energy function of the model consists of both 
a local and a nonlocal terms: 



H[<t>] = H L [<j>] + H NL [<fi] 



(4) 



According to <[3j for the case of zero external field, the local 
potential can be written in full generality: 

HM = \J rf 2 x[«(V0(x)) 2 - r0 2 (x) + ^ 4 (x)] (5) 

In order that the local potential assume the desired double- 
well structure, u and r are taken to be positive phenomeno- 
logical constants. We have also added a continuous approxi- 
mation for the exchange interaction in the form of an attractive 
square gradient term, which favors spatial homogeneity of the 
order parameter. 

Disregarding higher order interactions in (0, the nonlocal 
term modeling a repulsive dipolar interaction has the form: 



H NL [4>] 



^ J d 2 x y>x>(x)J(|x-x'|)0(x') (6) 



where J(|x — x'|) = l/|x — x'| 3 in two dimensions. In 
equations (0) and (O, k and 6 are positive phenomenological 
constants that describe, respectively, the range of the short- 
range exchange interaction and the strength of the long-range 
dipolar one. In all space integrations it is assumed a short 
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distance cutoff 2n/A, which in the numerical implementation 
appears naturally as the lattice constant a. 

Hence, in the limit of strong perpendicular anisotropy and 
strong damping, the relaxational dynamics of the system can 
be modeled by the following Langevin equations: 



dcj){x, t) 

Ft 



-r 



7?(x,t) 



(7) 



where T is a constant (the mobility) which sets the time scale, 
and rj(x, t) is a Gaussian thermal noise with (i](x, t)) = 
and {r](x,t)rj(x',t')) = 2 TTS(t - £') 5 2 {x - x'), as usual. 
Equation (0 can be explicitly written as: 



1 d(j>{x, i) 

f Ft 



«;V^(x,i) + rcf>(x,t) - u0 3 (x,t) (8) 
i /" d 2 x^(x',t)J(\x-x'\)+ V (x,t)/T 



It is convenient to express the above equations in dimen- 
sionless form by means of a transformation of variables. Fol- 
lowing Roland and Desai 20 , this transformation leads to the 
dimensionless form of equation dS): 



V 2 0(x,t)+0(x, t)-q) 3 (x,t) 
' f d 2 x'<f>(x',t)J{\x-x'\) + r)[x,t) (9) 



d^(x,t) 
dt 



with a short distance cutoff 2ir/A and thermal noise correla- 
tion (r)(x,t)r](x',t')) = 2TS(t - t')S 2 {x - x'). Now the 
parameter 6 stands for the relative strength between the two 
competing interactions. The last Langevin equation can be 
written in Fourier space as: 



d<f>(k,t) 
dt 



-A(k)q>(k, t) + [-0 3 (x, t) + V (x, t)] I (10) 



To be compatible with the subsequent numerical implemen- 
tation of this equation, the last two terms were not explicitly 
written in Fourier space. Here J means the k component of 
the corresponding Fourier transform. The function A(k) cor- 
responds to: 



A(k) 



1 + J(k)/S 



(11) 



which encodes all spatial information about the interactions. 
If this quantity has a negative minimum at a wave vector k m , 
selected by varying 5, the solution of equation ( fTOb is a mod- 
ulation in a single direction with periodicity given by k m . 

We have numerically solved the stochastic differential 
equation (fTOb discretizing space in a square lattice with mesh 
size a. Periodic boundary conditions were implemented us- 
ing the Ewald summation technique 22 in the long-range dipo- 
lar interaction. After spatial discretization this interaction is 
no longer isotropic for all spatial scales and it becomes grad- 
ually anisotropic as the wave vector comes close to it. At 
the relevant spatial scales for our simulations, J(k) is slightly 
anisotropic. It is important to note that symmetry properties 



of the different magnetic structures appearing at low temper- 
atures are affected by the square symmetry of the lattice. In a 
triangular lattice, for example, the phenomenology may be to 
some extent different^ 3 .. 

We found it advantageous to use a spectral method, since 
in the Fourier space form of the Langevin equation ( fTOb both 
spatial derivatives and the dipolar interactions acquire an al- 
gebraic form. 

The time derivative was approximated using a simple Euler 
scheme with a time step At. Taking an isotropic form of the 
discretized laplacian, the spatial derivatives were treated using 
a semi-implicit method, where the k 2 term in (fTOb is evaluated 
in the new time value. This treatment is standard to improve 
the stability of the algorithm 2 !. 

Therefore, discretizing equation (TToT > in such manner, 
i.e. through a first order semi-implicit spectral method, we ob- 
tain the following recurrence relation: 

0(k,t+At) = 1 - ^ fc2 { [1 + At- Ata 2 J(k)/6] <Xk, t) 
+ [-At<p\x 1 t)+T 1 {x 1 t)\l} (12) 

After discretization the noise term rj, in the way it appears 
in the last equation, is a random Gaussian number with ampli- 
tude (2T At/a 2 ) 1/2 . The dipolar interaction J(k) in CG) is 
the fast Fourier transform of the result of the Ewald summa- 
tion, evaluated at the beginning of the simulation. 

The computational advantage of updating the system in 
Fourier space is accomplished using an adaptative FFT 
algorithm 25 , where the main time consuming operations are 
transforming Fourier the field <fr and the (— (f> 3 + rj) term and 
then transforming back the new field value. 



III. RESULTS 

We performed simulations of the continuum dipolar model 
through equation ( fT2l . We analyzed the case 5 = 2\/5, where 
the ground state is a stripe modulated state with wave vector 
k m = 1.01447, which is close to ir/3. Simulations were per- 
formed for system sizes L = 192 and L = 384 for different 
temperatures. 

The existence of a nontrivial solution with wave vector k m 
of the Langevin equation ( fTOb at zero temperature depends on 
whether the term in (TTTb is negative at A: = k m . Before the 
variables transformation, this could be done varying the pa- 
rameter r. After the transformation, the existence of a solution 
is achieved only by tuning the mesh size a. We set a = 
The time integration was stable inside the range of tempera- 
tures of the simulations using a time step At = 0.5. 

To measure the orientational order of the stripes we first 
consider the director field: 



n(x) = 



V0(x) 



(13) 



which provides the local orientation of the stripe, when lo- 
calized on a domain wall between opposite magnetizations. 
An analogy between this quantity and the Frank director of 
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smectic liquid crystals has already been made on previous 
works 26 ' 27 ' 28 , and leads to the definition of a tensor order pa- 
rameter: 

Q a /3( X ) = «a( x )"-/3( x ) ~ 2 5a P ^ 

where a, (3 — 1,2 are the cartesian components. Similar 
to a nematic order parameter in a liquid crystal, an orienta- 
tional order parameter Q can be defined as the positive eigen- 
value of the spatial average of the above tensor order param- 
eter—. This value corresponds to the quantity cos 28, where 
n = (cos 9, sin (9) and 9 is the angle between the local direc- 
tor field and the mean orientation of the system. But, since 
the director field inside the stripe domains does not necessar- 
ily provide the local orientation of the stripe, to get a more 
precise value for the orientational order parameter we average 
Eq.(fl4b only over domain wall sites, namely: 

5a/S = ^XN«0(*) d5) 
x 

where the prime denotes the restricted sum with L' being the 
total number of domain wall sites. Now we can write explic- 
itly the orientational order parameter as: 

Q = y / Qn + Q?2 (16) 

One possible way to characterize the translational order of 
the stripes is through a staggered magnetization, defined as: 

m k = ^ sgn (0(x) ) sgn (cos(k • x))^ (17) 

where sng is the sign function and for k we use the ground 
state wave vector. 

The results we present here were obtained with the follow- 
ing procedure: the system is initialized in the ground state 
and heated with a heating rate dr- At each temperature of in- 
terest, the heating process is halted and the system is left to 
run a transient period of rit time steps before we start record- 
ing system configurations at each n m time steps, used later to 
measure the desired quantities. Typical values of rit were be- 
tween 5 x 10 4 to 4 x 10 6 in order to get as close as possible to 
equilibrium. To estimate both transient (n t ) and decorrelation 
(n m ) times we analyzed the behavior of the two-times correla- 
tion function for the different system sizes and temperatures. 
Typical n m values were in the order of 10 5 and the results are 
averages over 20 to 60 system configurations. 

In Fig. [T] we show the results for the translational and ori- 
entational order parameters. The data clearly shows a phase 
transition, where both kinds of order emerge at a finite tem- 
perature. In the L = 192 case we found that the two or- 
der parameters decay simultaneously to zero at T as 0.12, 
even though the staggered magnetization drops significantly 
at T = 0.11. In the region in between the translational order 
parameter seems to present strong finite size effects. As the 
temperature is lowered in the ordered phase, a stripe structure 
develops gradually by the annihilation of defects (see Fig. [3]). 
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FIG. 1 : Order parameters as a function of temperature for L — 192. 
The open symbols corresponds to the larger system size, L = 384. 
Thin vertical lines separate approximately temperature ranges of the 
different regimes. 

Although we were unable to study in more detail the finite 
size scaling near the transition, it seems improbable that an in- 
termediate, nematic-like phase, with only orientational order 
can be present in this model. Although translational order de- 
cays faster than the orientational one, our results point to the 
presence of a single phase transition from a disordered high 
temperature phase to a low temperature phase with both ori- 
entational and translational orders. This is one of the possible 
scenarios that emerge from a theoretical analysis of a similar 
model by Abanov et al. — . 

To take a closer look at the structural properties of the dif- 
ferent magnetic configurations, we calculated the static struc- 
ture factor: 

5(k)^(|0(k)| 2 ) (18) 
and the associated spatial correlation function, given by: 

k 

that can be quickly computed using a fast Fourier transform 25 . 
The relevant directions for the spatial correlation function 
are the directions parallel and perpendicular to the stripes, 
respectively denoted by C x and C y . This quantities de- 
scribe the translational order along the two relevant directions. 
We have also computed nematic (or orientational) correlation 
function a 27 ' 28 , since they encode information on the spatial de- 
cay of orientational order of the stripes. Nematic correlations 
are defined as: 

C nn (r) = ^^(TrQ(x + r)Q(x)) (20) 

X 
X 

+ Q 12 (x + r)Q 12 (x)> (21) 

where the function summed up in Eq. (|20| > is analogous to 
(cos [20(x + r) — 20(x)]). Examples of translational and ori- 
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FIG. 2: Correlation functions in the smectic-like region at T — 0.1. 
The upper figure shows the connected orientational correlation func- 
tion in the perpendicular direction for L — 384, plotted in log-log 
scale, before and after the smoothing procedure. The continuous line 
corresponds to a fitting by the function Ar~ OLy , with a y = 1.18 
and A a fitting parameter. In the parallel direction this function 
decays immediately to zero. The lower figure shows the paral- 
lel (inset) and perpendicular spatial correlations functions, both for 
L — 192. The parallel C x function is plotted in a log-log scale. The 
continuous lines corresponds to fittings by the functions r~" x and 
cos(kor)r~ u ' y with lu x = 0.33 and uo y — 0.19 and fco = 1.014. 
Similar values where found in the L — 384 case. 

entational correlation functions can be seen in Figs. [2] and [4] 
for two different temperatures. 

For low enough temperatures (T < 0.11), the nematic cor- 
relations are strongly affected by residual oscillations of the 
field 4>(x). In order to obtain accurate values for the orien- 
tational correlation function we found convenient to smooth 
the tensor order parameter ( TT4b following a smoothing proce- 
dure introduced in Ref. 28 . We smooth the fields Qn(x) and 
Qi2(x) using for each the iterative process: 

/(n+l)(x) = i/(„)(x) + ^ Z(«)( X ') < 22 ) 

x'eAW 

where f/ n -\ is one of the fields after n iterations, and NN 
means the four nearest neighbors of x on the square lattice. 
We found that three iterations are enough to get sensible re- 
sults (see Fig|2). Similar smoothing procedures for orienta- 
tional correlation functions were used in simulations of melt- 
ing of two-dimensional solid system s 30 ' 31 . For clarity, in this 
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at (a) T = 0.1 and (b) T = 0.104. In (a) the encircled region 
exhibits two pairs of dislocations; above there is a "passage" due 
to two dislocations inserted in the same stripe, and below a pair of 
dislocations separated by one stripe. In (b) a pair separated by four 
stripes and a pair formed by a larger Burgers vector dislocation and 
a double dislocation. 



work we always show the connected orientational correlation 
function, where the mean square orientational order parame- 
ter is substracted from Eq.d2lY In this way, C nn accounts for 
direction fluctuations only. 

At sufficiently low temperatures (T < 0.03) spatial corre- 
lations decay immediately to a constant and there is transla- 
tional and orientational long-range order (see Fig. Q] for the 
values of the order parameters in the corresponding tempera- 
ture regimes). 

At slightly higher temperature and in a wide range (0.03 < 
T < 0.11) we found evidence of a low temperature smectic- 
like regime, similar to the smectic crystal phase predicted by 
Abanov et alJi and observed experimentally by Portmann et 
alA In this phase, undulation fluctuations (meandering exci- 
tations) alone are sufficient to cause algebraic decay of trans- 
lational correlations at low temperatures for T > 0.03. But 
at T = 0.09 bound pairs of dislocations appear and become 
more common as temperature increases. Long-range orienta- 
tional order persists over a wider range of temperatures, but 
above T = 0.065 orientational order also starts to decay alge- 
braically. 

An example of the algebraic behavior of the correlation 
functions in this region is shown in Fig. [2] for T = 0.1. It 
is important to note that the exponent of the translational al- 
gebraic decay in this regime increases with temperature and 
is different in the perpendicular and parallel directions - ex- 
ponents values range from u> x — 0.034 and to y = 0.036 for 
T = 0.035 to u> x = 0.365 and uo y = 0.212 for T = 0.104. As 
for the orientational order we observe a small temperature de- 
pendence of the algebraic decay exponent in the perpendicular 
direction, with exponents lying in the range 1.18 < a < 1.43. 
In the parallel direction, the orientational correlation function 
is a constant, indicating long-range order in this direction. 

In Fig. |3] we illustrate some typical configurations of this 
phase. We see that it is characterized by undulation excitations 
and a finite density of dislocation pairs - some different types 
are shown encircled. 
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FIG. 4: Correlation functions in the nematic-like region at T = 
0.114 and L = 384. The upper figure shows the connected ori- 
entational correlation function in the perpendicular and parallel di- 
rections in a log-linear scale. The continuous lines corresponds to 
fittings by the function Aexp(—r/\), with X y — 13.77 and X x = 
6.42. The lower figure shows the parallel (inset) and perpendicu- 
lar spatial correlations functions. The parallel C x function is plotted 
in a log-linear scale. The continuous lines correspond to fittings by 
the functions Aexp(—r /£ x ) and cos(fcof) exp(— r/£ y )r~" y with 
£ x = 10.17 and £ y = 60.0 and cu y = 0.26 and k = 0.98. 



In a narrow range of temperatures, 0.11 < T < 0.118, were 
the orientational order parameter is still high but the staggered 
magnetization drops considerably fast, we found a change of 
behavior in the correlation functions. In the parallel direc- 
tion, translational order is decorrelated exponentially rapidly 
to zero, indicating absence of translational order of the stripes. 
In the perpendicular direction, where order is more robust, 
there is an intermediate kind of behavior where the transla- 
tional correlation function is better fitted by a product of a 
power law and an exponential. An example of this behavior 
is shown in Fig. [4] where one can also see that the connected 
orientational correlation functions now decay exponentially in 
both directions but with different correlation lengths. Orien- 
tational domains in this regime are larger in the perpendicular 
direction. 

Analyzing visually the configurations of this region we 
observe many excitations disrupting orientational order. A 
typical configuration is shown in Fig. [5] We see that this 
regime presents large Burgers vector dislocations that can 
be regarded as a tightly-bound pair of oppositely charged 
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(a) 




FIG. 5: Typical configuration in the nematic-like region at (a) T = 
0.114 and (b) T = 0.116 for L = 384. In (a) some defects are 
put in evidence. In small circles, two pairs of dislocations inside 
smectic-like domains; in the larger circle a disclination dipole; shown 
colored a series of dislocations inside a single stripe and in an ellipse 
a domain of perpendicular orientation. 



disclinations . There are also what may be called a dislo- 
cation cascade, a series of bifurcations within a single stripe 
(colored in Fig. [5^). Small domains of perpendicular ori- 
entation (encircled in Fig. ((5^)) are present as well, and be- 
come more common and larger as temperature increases (see 
Fig. \5fa). Surrounded by topological defects there are do- 
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FIG. 6: Orientational correlation length (full symbols) and suscepti- 
bility (open symbols) as a function of temperature in the isotropic 
phase. The continuous curves show power-law fittings, namely: 
£nn(T) ~ (T — 0.1172)" 1 - 19 and X nn(T) ~ (T — 0.1161)- ' 61 . 



mains of locally smectic-like arranged stripes, where order is 
decorrelated mainly by meandering excitations and all kinds 
of pairs of dislocations (some are encircled in Fig.[5k). 

In order to estimate the transition temperature T c between 
the orientationally ordered and isotropic phases, we have mea- 
sured the orientational correlation length and susceptibility of 
Q in the isotropic phase. The results for L — 192 are shown 
in Fig. [6], together with power-law fits. We found that T c lies 
close to T c ps 0.117 and that divergences are well fitted by 
power-law forms, at variance with the exponential divergence 
that one would expect in a defect-mediated transition accord- 
ing to the KTHNY theory 26,30 . The power law behavior seems 
to be in agreement with a recent result in which the nematic 
transition is predicted to be second order 18 . 

Finally, we discuss the evolution of the structure factor, 
eq.([T8l>, with temperature. In Fig. Q we show four charac- 
teristic examples around the transition. From the first to the 
second plot it can be seen that the sharp peaks characteriz- 
ing stripe order are replaced by nematic-like peaks, spreaded 
along the ring |k| ~ k m due to angle fluctuations. The tran- 
sition to the disordered phase seems to be through an increase 
of domains of perpendicular orientation of the stripes. In the 
structure factor, this is reflected as the growth of two symmet- 
ric peaks around k i=s k m in the perpendicular direction. Im- 
mediately after the transition, at T = 0.118, the four symmet- 
rical peaks are equivalent, as can be seen in the third plot. As 
the temperature is further increased, angle fluctuations around 
this two preferential directions smear out the four peaks and 
the system becomes almost isotropic and the spectral weight 
of the structure factor lies on a ring with a weakly tetragonal 
shape, as shown in Fig. [7] A typical configuration illustrating 
the weakly tetragonal symmetry of the disordered phase just 
above the transition is shown in Fig. [8] 

Not far away from the transition, the isotropic phase still 
presents the exponential-algebraic behavior of positions cor- 
relations as in the nematic-like regime, indicating that it lo- 
cally resembles the low temperature phase. The anisotropic 
feature of the orientational correlation function dissapears in 




FIG. 7: Structure factor at T = 0.1, 0.114, 0.118 and 0.14 for L 
384. 



this phase, where C nn decays exponentially to zero. 



IV. CONCLUSIONS 

We have made Langevin simulations at finite temperature of 
a model for ultrathin ferromagnetic films with perpendicular 
anisotropy. These systems have a stripe ground state due to 
competition between exchange and dipolar interactions. We 
have found a phase transition from a high temperature dis- 
ordered phase to a low temperature phase with both trans- 
lational and orientational orders. Below the transition point 
isotropy is spontaneously broken, and a smectic-like magnetic 
structure develops. Both kinds of ordering seem to take place 
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FIG. 8: The left configuration shows the isotropic phase at T = 0.12 
for L = 192, presenting a target defect close to its center. Small 
target defects are ubiquitous in this phase. At top right we show a 
large target defect in a portion of a L = 480 configuration in an 
out-of-equilibrium disordering process at T — 0.13. In the middle 
a portion of a L = 384 configuration at T — 0.112 exhibiting a 
large dislocation pair made of disclination pairs. Below, pairs of 
dislocations at T = 0.104. 

at the same T c , in agreement with some previous theoretical 
predictions. The absence of an intermediate, purely orien- 
tational, phase is probably due to the isotropic nature of the 
model Hamiltonian. Other terms, which explicitly break ori- 
entational order, may be necessary in order to have an inter- 
mediate nematic phase. This is what emerges from a recent 
analysis of the nematic transition in this same model, where 
an extra term was included and shown to be responsible for 
the existence of an intermediate nematic phase 18 . The nature 
of these terms depend on microscopic interactions, namely, on 
the modelling of the anisotropy contribution. In this respect, 
more input from experiments is fundamental. 

We have found that the numerical results show complex 
magnetic structures, with translational and orientational cor- 
relations decaying algebraically at low temperatures. Transla- 
tional correlations have different exponents for the longitudi- 
nal and perpendicular components, relative to the ground state 
stripe orientation. The exponents are also temperature depen- 
dent. In the case of orientational correlations, we observe a 
weak temperature dependence in the perpendicular direction, 
and a saturation to a constant value in the longitudinal direc- 
tion. Qualitatively one can say that we observe translational 
quasi-long-range order and orientational long range order, at 
low temperatures. 

As temperature grows it is possible to observe a prolifer- 
ation of topological defects, with structures similar to those 
observed in real systems. In a narrow region around the tran- 
sition temperature, both translational and orientational corrre- 



lations begin to decay exponentially. A power law fit of the 
data of the orientational correlation length and the associated 
susceptibility in the high temperature region allows us to es- 
timate the critical temperature. This power law dependence 
of the orientational correlation length does not agree with the 
well known KTHNY theory of two dimensional melting. We 
do not know to what extent the predictions from this theory 
should be applicable to magnetic two dimensional systems 
like this. Our results are in agreement with a possible sce- 
nario for a smectic magnetic system put forward by Abanov 
and coworkers more than 10 years agoii, and with a recent 
result on the second order nature of the nematic transition in 
the same model studied in this work supplemented by another 
term in the Hamiltonian which explicitly breaks orientational 
order 18 . 

We found that the weak anisotropy of the dipolar interaction 
due to spatial discretization and the short width of the stripes 
(of three grid points) has led to a pinning effect that favored 
the stripes orientation to be preferentially on the two Cartesian 
directions. But looking at the structure factor around the tran- 
sition, we found that the fluctuations responsible for disrupt- 
ing orientational order are not restricted to the two Cartesian 
directions and the isotropic nature of the model still manifests. 

More quantitative analytic predictions are clearly needed 
in order to assess the quality and limitations of our results. 
Also, it would be extremely interesting to do systematic ex- 
perimental measurements of structure factors and correlations 
as a function of temperature, like the ones done by Pescia and 
coworkers 2 ' 10 ' 3 2 . This would allow, between other things, 
to elucidate the experimental conditions under which a ne- 
matic intermediate phase can be present in a magnetic sys- 
tem. The similarity between our simulations and experimental 
images of perpendicularly magnetized fee Fe films grown on 
Cu(100) 2 ' 10 ' 32 is striking. For example, target defects first ob- 
served by Vaterlaus et al. 10 were observed in our simulations 
and are ubiquitous in the disordered phase (see Fig. [8). 

Finally, the dynamical behavior of these systems also de- 
serve to be studied, because the presence of frustration and 
the emergence of many kinds of topological defects lead to a 
complex dynamics, with pinning of magnetic structures dur- 
ing long time scales, before the final, asymptotic equilibrium 
state sets ini2. 
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